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Abstract 


The  MOS  transistor  physical  model  as  described  in  [3]  is  presented  here  as  a  network  model.  The  goal  is 
to  obtain  an  accurate  model,  suitable  for  simulation,  free  from  certain  problems  reported  in  the  literature 
[13],  and  conceptually  as  simple  as  possible.  To  achieve  this  goal  the  original  model  had  to  be  extended 
and' modified.  The  paper  presents  the  derivation  of  the  network  model  from  physical  equations,  including 
the  corrections  which  are  required  for  simulation  and  which  compensate  for  simplifications  introduced  in 
the  original  physical  model. 

Our  intrinsic  MOS  model  consists  of  three  nonlinear  voltage-controlled  capacitors  and  a  dependent  current 
source.  The  charges  of  the  capacitors  and  the  current  of  the  current  source  are  functions  of  the  voltages 
y  Vbs,  and  Vds.  The  complete  model  consists  of  the  intrinsic  model  plus  the  parasitics.  The  apparent 
simplicity  of  the  model  is  a  result  of  hiding  information  in  the  characteristics  of  the  nonlinear  components. 
The  resulted  network  model  has  been  checked  by  simulation  and  analysis.  It  is  shown  that  the  network 
model  is  suitable  for  simulation:  It  is  defined  for  any  value  of  the  voltages;  the  functions  involved  are 
continuous  and  satisfy  Lipschitz  conditions  with  no  jumps  at  region  boundaries;  Derivatives  have  been 
computed  symbolically  and  are  available  for  use  by  the  Newton-Raphson  method.  The  model’s  functions 
can  be  measured  from  the  terminals.  It  is  also  shown  that  small  channel  effects  can  be  included  in  the 
model.  Higher  frequency  effects  can  be  modeled  by  using  a  network  consisting  of  several  sections  of  the 
basic  lumped  model. 

Future  plans  include  a  detailed  comparison  of  the  network  model  with  models  such  as  SPICE  level  3  and 
a  comparison  of  the  multi-section  higher  frequency  model  with  experiments. 
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1  Introduction 


Our  interest  in  modeling  MOS  transistors  started  from  our  work  on  the 
Supercomputer  Toolkit  [1,  2].  The  Supercomputer  Toolkit  is  a  family  of 
hardware  modules  (processors,  memory,  interconnect,  and  input-output  de¬ 
vices)  and  a  collection  of  software  modules  (compilers,  simulators,  scientific 
libraries,  and  high-level  front  ends)  from  which  high-performance,  special- 
purpose  computing  systems  can  be  easily  configured  and  programmed.  The 
main  advantage  of  the  Supercomputer  Toolkit  is  that  for  numerically  inten¬ 
sive  computations  it  provides  supercomputer  class  performance  at  a  fraction 
of  a  supercomputer  cost.  One  such  application  is  the  simulation  and  opti¬ 
mization  of  high  speed  circuits  with  respect  to  physical  parameters. 

For  optimization  of  Bi-CMOS  circuits  we  needed  a  model  of  the  MOS 
transistor.  Our  colleagues  (see  acknowledgment)  recommended  a  long  list 
of  references,  including  [3,  4,  5,  6].  Tsividis  [3]  is  the  most  comprehensive 
book  of  its  kind.  The  simplest  approach  would  have  been  to  use  the  SPICE 
level  3  model  [7].  However,  several  shortcomings  of  this  model  have  been 
reported  [13,  12],  all  of  them  involving  the  simulation  of  analog  circuits  of 
the  type  we  intended  to  use.  With  the  Toolkit  in  hand,  computation  time 
did  not  seem  the  hindering  factor  for  the  problems  we  intend  to  attach,  so 
we  decided  to  construct  the  most  accurate  model  possible.  Thus,  the  above 
references,  in  particular  [3],  became  the  basis  for  this  work. 

The  type  of  transistor  model  we  wanted  was  a  network  model ,  i.e.,  one 
that  models  the  transistor  as  a  network  of  more  primitive  devices.  The  re¬ 
sulting  network  can  be  use  as  a  component,  i.e.,  interconnected  with  other 
components  to  form  larger  networks.  The  appropriate  transistor  model  de¬ 
pends  on  the  nature  of  the  required  solution.  In  particular,  for  transient 
simulation  the  model,  once  embedded  in  a  network  satisfying  certain  con¬ 
ditions,  should  lead  to  an  appropriate  set  of  ordinary  differential  equations 
which  can  be  integrated  numerically  to  yield  a  solution  (see,  for  example, 
110])- 

The  material  we  found  in  the  above  references  was  not  quite  in  the  form 
we  wanted.  There  was  no  large  signal  network  model  and  only  the  small 
signal  model  could  be  described  as  a  network  model.  The  number  of  details 
was  very  large.  The  emphasis  was  on  deriving  the  model  from  physical 
considerations  and  we  missed  a  structural,  hierarchical  description  of  the 
results  where  the  details,  if  needed,  could  be  looked  up  later  and  filled  in.  The 
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physical  equations  were  derived  as  if  the  transistor  had  infinite  dimensions. 
The  lumped  models  derived  from  these  equations  are  approximations  to  the 
behavior  of  a  finite  size  transistor.  Yet,  we  could  not  find  a  mention  of  what 
properties  could  or  could  not  be  measured  from  the  terminals.  Thus,  our 
goal  was  to  start  from  physical  considerations  (using  the  above  references) 
and  find  a  “simple”  large  signal  network  model  whose  components  can  be 
measured  from  the  terminals. 

In  the  next  section  we  derive  the  model  for  the  usual  operating  region 
of  the  MOS  transistor.  The  result  is  presented  in  figure  2.  Section  3  gives 
our  reasons  for  extending  the  model  to  other  operating  regions  and  then 
presents  such  extension.  Section  4  presents  the  results  graphically.  Section 
5  addresses  the  issue  of  measuring  the  characteristics  externally.  Section  6 
discusses  some  properties  of  the  model,  in  particular  those  related  to  the 
existence  and  uniqueness  of  the  response  of  networks  constructed  with  our 
transistor  models  as  components.  The  last  section  presents  the  summary  and 
the  conclusions. 

Appendix  A  presents  the  model  equations  in  the  “usual”  operating  region. 
Appendix  B  discusses  the  derivation  of  the  incremental  capacitance  from  the 
charge  functions  using  symbolic  computing  techniques. 

2  Region  I  -  The  usual  operating  region 

This  section  presents  the  model  in  the  usual  operating  region.  The  notations 
follows  [3].  The  treatment  is  brief  and  justification  can  be  found  in  the  above 
reference.  Figure  1  presents  the  transistor  including  the  definitions  of  the 
voltages  and  currents  at  its  terminals. 

Modeling  starts  by  considering  the  behavior  of  the  transistor  in  several 
regions  defined  by  the  terminal  voltages  Vgs ,  14,  and  V*.  Region  I  is  defined 
by  the  following  relations  among  the  voltages: 

14,  >o 

14  5:  0. 

This  region  is,  in  turn,  divided  into  four  subregions: 

Accumulation:  Vgs  <  Vfb  +  14,  where  Vfb  is  the  flat  band  voltage; 

Cutoff:  Vfb  +  14s  <  Vgs  <  Vth,  where  Vtb  is  the  cut-off  voltage; 
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Figure  1:  The  MOS  transistor,  it’s  terminal  voltages  and  currents,  and  a 
differential  cross  section. 

Saturation:  Vgs  >  Vth  and  Vds  >  Vdsat ,  where  Vdsat  is  the  saturation  voltage; 
Linear:  Vgs  >  Vth  and  Vds  <  Vdsat. 

In  the  above,  Vth  is  the  threshold  voltage,  a  function  of  Vj>4, 

Vth  =  Vto  +  7  *  {\J fys  —  Vbs  —  \f$s), 
and  Vdsat  i  the  saturation  voltage,  is 

T,  Vgs  -  Vth 

Vdsat  —  i  i  r 

l  +  o 


2^1  +  4>s-  U,  ’ 

where  Vto,  7  and  are  physical  constants  (i.e.,  independent  of  the  volt¬ 
ages)  which  are  defined  in  Appendix  A. 

To  compute  the  current  and  the  charges  one  considers  a  differential  sub¬ 
section  of  the  transistor  (see  Figure  1).  Assuming  the  fields  in  the  transistor 
are  the  same  as  if  the  transistor  were  infinite  in  the  x  and  y  direction,  from 
physical  consideration,  one  calculates  the  current  through  the  subsection  and 
the  charges  in  it  as  functions  of  certain  physical  parameters  and  of  the  field 
intensity  in  the  subsection.  This  field  intensity  can  be  expressed  in  terms  of 
the  voltages  Vgs,  Vbs ,  and  Vds. 
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The  current,  I, is .  is  the  crrent  flowing  through  the  differential  cross  sec¬ 
tion,  in  the  (minus)  x  direction  due  to  both  drift  and  diffusion.  It  is  assumed 
that  only  a  negligible  amount  of  this  current  enters  the  bulk  and,  therefore, 
the  current  entering  the  cross-section  on  the  left  exits  on  the  right.  From 
that  also  follows  that  its  value  can  be  obtained  by  integration  along  the  the 
subsection  face  (integration  along  the  0  axis  and  the  y  axis;  the  distribution 
along  the  y  axis  being  uniform  as  a  result  of  the  infinite  dimension  assump¬ 
tion).  The  integration  results  in  the  following  expression  for  the  linear  and 
saturation  regions: 


Ids  —  2 L^CoxVdsatiVgs  ~  ft/i)(l  —  <*2)(1  +  AVrfa)  (1) 

where  W  and  L  are  the  width  and  length  of  the  transistor,  Cox  is  the  gate 
oxide  capacitance  per  unit  area,  a  is  zero  in  the  saturation  region  and  equal 
to 

x  Vu 

Vdsat 

in  the  linear  region.  The  last  term,  (1  +  A 14,)  is  (alas!)  a  correction  term 
introduced  experimentally  (it  does  not  follow  from  the  physical  considerations 
in  the  infinite  dimension  x  ). 

As  we  move  from  saturation  to  cutoff  the  transistor  passes  through  a 
medium  and  then  a  weak-inversion  zone;  Ids  goes  down  sharply  and  eventu¬ 
ally  diminishes  to  zero: 


Ids  =  Io  exp 


Vga  -  Vth 

3  kT  ’ 


2  q 


(2) 


where 


,  _  W  _  3  ,kT  2 

Io  -  L  yCox^(  q  )  e  . 


To  make  Ids  a  continuous  function  of  Vgs,  I0  has  to  be  added  to  the  value 
of  Ids  in  the  saturation  and  linear  regions  (see  appendix  A  for  the  complete 
equation). 

The  charges  provide  a  more  intriguing  story.  In  the  infinitesimal  section 
dx  above,  four  charges  (per  unit  length)  are  defined: 


•  qg  -  the  charge  on  the  metal  gate; 
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•  qi  -  the  inversion  layer  charge; 

•  qb  ~  the  charge  in  the  bulk; 

•  qo  -  the  interface  charge,  a  fixed  quantity  independent  of  the  fields  in 
the  section. 

These  charges  satisfy 

qg  +  qi  +  qb  +  qo  —  0. 

Each  of  these  charges  can  be  expressed  as  a  function  of  the  position  x 
and  the  voltages  Vgs ,  Vbs  and  Vds •  The  charges  can  also  be  integrated  over 
x  to  yield  the  total  gate  charge,  Qg ,  the  total  inversion  (layer)  charge,  Q/, 
the  total  bulk  charge,  Qb ,  and  the  total  interface  charge,  Q0-  The  sum  of  the 
total  charges  is  still  zero,  and  so  is  the  sum  of  the  charges’  derivatives  with 
respect  to  time. 

The  above  relations  for  the  charges  have  been  derived  under  DC  con¬ 
ditions.  It  is  assumed  (the  pseudo-static  assumption)  that  as  the  voltages 
change,  but  “not  too  fast,”  the  charges  change  according  to  the  statically  de¬ 
rived  relations.  (As  usual,  a  quantitative  value  for  “not  too  fast”  is  not  given 
in  the  literature  and  frequency  bounds  on  the  validity  of  the  pseudo-static 
model  have  to  be  obtained  experimentally.) 

The  next  step  is  conceptually  delicate:  We  would  like  to  partition  Qj 
into  two  charges  Qd  and  such  that  Qi  =  Qs  +  Qd',  and,  to  associate  each 
charge  with  the  corresponding  terminal  —  Qs  with  the  source,  Qd  with  the 
drain,  Qg  with  the  gate  and  Qb  with  the  bulk.  Moreover,  we  claim  that,  given 
that  the  proper  partition  is  done,  the  time  derivative  of  each  of  these  charges 
is  (a  component  of)  the  current  in  the  corresponding  terminal,  i.e.  =  igi 
etc.  Except  for  the  case  of  Qg  the  support  of  this  step  is  not  trivial;  it  is 
given  in  detail,  together  with  a  proper  partition,  in  chapter  7  of  [3]  and  its 
references.  One  of  the  arguments  given  there  is  that  the  results  agree  with 
experiment. 

Notice  that  since 

Qg  d-  Qi  T  Qb  T  Qo  —  0, 

and  ftQ0  =  0, 

d  ^  <:/  f/  d 

JtQg  +  jtQs  +  jtQa  +  jtQb  =  0.  (3) 
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In  words,  the  above  equation  means  that  the  terminal  currents  due  to  changes 
in  the  charges  satisfy  Kirchoff’s  current  law. 

Thus,  we  may  model  the  effect  of  the  charges  and  Ids  using  three  voltage- 
controlled  non-linear  capacitors1  and  a  dependent  current  source  (see  Figure 
2).  The  characteristics  of  Cga,  for  example,  is  given  by  the  function  specifying 
Qg  as  a  function  of  Vgs,  Vb3,  and  Vds ■  This  means  that  this  capacitor  is  a 
four-terminal  component,  needing  a  reference  and  three  terminals  to  specify 
its  charge  (see  Figure  3.  The  time  derivative  of  the  charge  on  Cgs ,  the  charge 
Qg ,  is  a  current  entering  g  and  leaving  through  s.  Thus,  changes  in,  say,  Vds , 
may  change  Qg  and  cause  a  gs  current.  Cb3  and  Cds  are  similarly  described 
-  each  with  its  own  characteristics.  The  functions  are  given  in  Appendix  A. 
The  characteristics  of  the  current  source  is  described  by  equations  1  and  2. 

The  choice  of  the  non-linear  capacitor,  rather  than  the  commonly  chosen 
incremental  capacitor,  is  a  major  one.  It  provides  a  simple  link  with  the 
charge  equations  and,  at  the  same  time,  simplifies  the  model  conceptually. 
The  incremental  capacitors  are  the  partial  derivatives  of  the  characteristics 
at  a  given  operating  point. 

Note  that  we  can  define  a  capacitor  for  each  of  the  terminal  charges. 
Any  three  of  the  four  can  be  chosen  for  the  model.  The  network,  or  rather 
Kirchoff ’s  law,  guarantees  that  the  sum  of  the  terminal  currents  is  zero.  From 
equation  3  follows  that  the  omitted  charge  is  equal  to  the  minus  of  the  sum 
of  the  other  three  charges  (plus  a  constant  which  we  are  at  liberty  to  take  as 
Qo).  Note  also  that  in  Figure  2  the  source  s  was  chosen  as  a  reference  node 
and  the  charge,  say,  Qg  was  associated  with  the  capacitor  Cgs.  This  choice 
yields  characteristics  whose  expressions  corresponds  directly  to  the  charge 
equations  in  the  literature  [3].  If,  however,  the  base  or  gate  is  chosen  as  the 
reference  node,  symmetry  implies  that  the  characteristics  of  Csb  is  equal  to 
the  one  of  C,ib  ■ 

The  model  of  Figure  2  is  called  the  the  intrinsic  model.  To  this  model 
we  have  to  add  the  so-called  “parasitics” .  These  are  two  linear  capacitors 
representing  the  overlap  between  the  gate  metallic  layer  on  both  the  source 
and  the  drain,  respectively,  and  two  nonlinear  capacitors  representing  the 
effect  of  the  pn  junctions  between  the  source  to  bulk  and  drain  to  bulk.  A 

1 A  capacitor  is  called  voltage-controlled  if  its  charge  is  given  as  a  function  of  its  voltage, 
i.e.,  q  =  f(v);  The  function  is  called  the  characteristics  of  the  capacitor.  In  our  case,  v 
is  a  vector  whose  components  are  voltages.  The  extension  of  the  terminology  is  straight 
forward. 
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Cgs 


Figure  2:  The  intrinsic  model 


Figure  3:  A  four-terminal  capacitor;  Q  is  a  function  of  (Vi ,  V2,  T3);  The 
current  i  =  ^Q. 
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Figure  4:  The  complete  model;  The  intrinsic  model  plus  parasitics. 


better  model  for  the  pn  junction  is  that  of  a  non-linear  capacitor  in  parallel 
with  a  non-linear  resistor  which  represents  a  diode;  In  region  I,  however,  the 
diode  is  reversed  biased,  so  the  resistor  can  be  neglected. 

Figure  4  describes  the  complete  model.  The  functions  describing  the 
characteristics  of  the  components  appear  in  Appendix  A. 

Discussion:  The  above  contains  an  outline  of  the  modeling  process. 
Most  of  the  details  are  tucked  away  in  the  functions  that  describe  the  com¬ 
ponent  characteristics.  We  would  like  to  consider  these  functions  from  the 
point  of  view  of  the  model  use,  i.e.,  circuit  simulation.  We  expected  the  func¬ 
tions  to  satisfy  Lipschitz  conditions  in  terms  of  the  independent  variables  2 
which,  here,  are  the  voltages.  Extracting  the  functions  from  literature  we 
came  across  several  surprises.  The  literature  provided  algebraic  expressions 
for  each  region  (i.e.,  the  saturation  region,  cutoff,  etc.)  but  the  functions  were 
not  continuous  at  the  boundaries.  This  we  discovered  by  displaying  cuts  of 
the  functions  at  various  values.  As  a  result,  the  functions  in  Appendix  A 
contain  corrections  that  we  introduced  heuristically. 

In  circuit  simulation,  the  Newton-Raphson  method  is  used  to  solve  the 
algebraic  equations  that  result  from  implicit  integration  methods.  Newton- 


Raphson  requires  the  partial 
voltages;  E.g.,  *,Vd^. 


derivatives  of  the  functions  with  respect  to  the 
Deriving  these  twelve  functions  was  not  simple. 


The  expressions  in  each  region  are  complicated  enough  so  we  just  could  not 


2  In  order  to  have  a  solution  defined  by  a  set  of  first  order  differential  equations  the 
network  has  to  satisfy  certain  conditions.  One  of  these  conditions  is  that  the  functions 
satisfy  the  Lipschitz  conditions  in  term  of  the  independent  variables. 


9 


Figure  5:  The  five  regions.  Lines  (except  the  V(i.?  coordinate)  are  boundaries 
between  regions. 

get  the  evaluation  right  when  we  did  the  differentiation  by  hand.  We  found 
that  we  had  a  wrong  function  on  hand  by  testing.  Every  solution,  i.e.,  pro¬ 
posed  derivative  function,  was  compared  with  the  numerically  approximated 
value;  when  the  discrepancy  was  too  large  it  indicated  a  mistake. 

Thus,  we  turned  to  symbolic  algebra  to  evaluate  the  derivatives.  The 
straightforward  derivation,  however,  produced  enormous  expressions  which 
we  could  not  simplify  well.  We  circumvented  this  difficulty  by  writing  a 
symbolic  program  for  evaluating  each  derivative.  This  technique  is  described 
in  Appendix  B  by  example.  Figure  18  depicts  the  partial  derivative  of  Qg 
with  respect  to  Vgs. 

It  is  worthwhile  noting  that  these  derivatives  have  circuit  meaning.  They 
are  the  incremental  capacitance  and  trans-capacitance  that  appear  in  the 
small  signal  model.  There  are  nine  of  those  —  the  partial  derivatives  of  Qg, 
Qbi  and  Qd  each  with  respect  to  each  of  the  voltages  Vgs,  Vj,*,  and  Vds. 

On  inspection  one  sees  that  these  derivatives  are  not  continuous.  The  dis¬ 
continuities  occur  along  the  subregions’  boundaries.  In  all  cases,  the  deriva¬ 
tive  normal  to  the  boundary  has  the  same  sign  on  both  side  of  the  boundary. 
This  fact  is  aesthetically  displeasing.  It  does  not,  however,  effect  the  exis¬ 
tence  and  uniqueness  of  the  transient  solution.  In  our  circuit  examples  it 
did  not  affect  the  convergence  of  the  Newton-Raphson  iteration.  It  does, 
however,  indicate  a  need  to  make  the  small  signal  model  more  accurate  at 
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Figure  6:  A  picturial  illustration  of  evaluating  the  functions  Qd  and  Qs  in 
each  region.  In  region  I,  Qd  and  Qs  are  evaluated  from  Vgs,  Vbs  and  Vds.  In 
region  II  the  V’s  are  evaluated  first.  This  amounts  to  translating  Vgs,  \ fbs 
and  Vds  of  region  II  to  a  point  in  region  I  illustrated  by  the  curved  arrow 
from  region  II  to  region  I.  Qd  and  Qs  are  evaluted  at  that  region  I  point;  The 
value  of  Qd  in  region  II  is  the  value  just  obtained  for  Qs  in  region  I.  In  other 
regions  the  functions  are  evaluated  similarly. 
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region  boundaries. 


3  The  Five  Regions  Model 

We  have  coded  the  Region  I  model  and  run  some  simple  examples.  Again, 
to  our  surprise,  the  simulation  trajectory  crept  out  of  Region  I  into  regions 
where  the  model  was  not  valid.  A  check  of  a  MOS  circuit  (An  inverter 
with  active  /^channel  load,  [11]  page  391)  showed  that  indeed  the  actual 
circuit  leaves  region  I.  The  phenomena  is  not  a  feature  of  the  model  or  of  the 
simulation  program  but  appears  in  actual  circuits.  Therefore,  we  constructed 
a  model  that  is  valid  for  all  values  of  Vgs,  Ha,  and  Vds-  We  refer  to  this  model 
as  the  five  regions  model.  Figure  5  gives  the  five  regions  that  partition  the 
(Vbs,  Vds)  space.  For  each  region  we  define  expressions  that  are  valid  for  any 
value  of  Vgs. 


3.1  Region  I 


We  return  to  region  I  in  order  to  extend  this  region  somewhat.  Considering 
the  functions  in  region  I  we  find  that  they  behave  properly  for  Vbs  >  0  except 
that  some  of  the  derivatives  become  infinite  (blow  up)  for  Vj,,  =  0,.  Thus, 
region  I  can  be  extended  to 

Vd.  >  0 

Vbs  ^  ^regions— border  , 

where  (fregions-border  <  4>s-  In  OUr  simulations  we  took  fregions-border  =0.50,. 

The  point  Vds  =  0,  Vgs  =  Vth,  i.e.,  Vdsat  =  0,  is  rather  delicate.  The  various 
functions  are  continuous  but  the  numerical  evaluation  is  apt  to  “blow  up”  as 
a  result  of  devision  by  zero  if  care  is  not  taken.  In  our  own  representation  of 
the  function  the  phenomina  occurs  in  the  calculation  of  alpha  as 

Vu 


a  =  1 
0 


Vdsat 


a 


for  Vds  <  Vdsat 
for  Vds  ^  Vdsat 


caused  a  division  by  zero  for  Vds  zero, 
eliminated  the  numerical  trouble: 

Vds 


A  slight  change  in  the  conditions 


a  =  1 


Vdsat 


foT  (Ids  ^  Vdsat) dTldi^Vdsat  ^  0) 
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a  =  0 


otherwise. 


The  fact  that  this  change  helped  depends  on  the  numerical  properties 
of  our  system,  in  particular,  the  way  the  system  rounds  a  small  number  to 
zero.  In  general,  note  that  at  the  above  point  functions  such  as  Qg  have  a 
double  zero  at  the  numerator  and  a  single  zero  at  the  denominator  and  the 
numerical  result  might  depend  on  the  order  of  computation. 

3.2  Region  II 

The  region  is  defined  by  the  symmetry  of  the  MOS  transistor.  Clearly,  the 
source  and  drain  are  symmetric  and  for  appropriate  voltages  the  source  acts 
as  drain  and  the  drain  acts  as  source.  Region  II  is  the  region  where  we  would 
be  in  region  I  if  we  interchanged  the  source  and  the  drain.  Let  us  name 
this  replacement  of  source  by  drain  “the  exchange  operation.”  Denote  the 
voltages  before  the  exchange  by  V,  and  the  voltages  after  the  exchange  by 
V .  The  relation  between  the  the  V’s  and  the  V’s  can  be  simply  concluded 
from  Figure  1: 

Vgs  —  Vgs  Vds  , 

Vbs  —  Vbs  Vds, 

Vds  =  -Vu. 

It  is  easy  to  verify  that  these  equations  map  the  following  region  into 
region  I: 


{ (Lgs?  14,  Vds )  | (14  ^  ^regions— border  >  Vds  V  Vbs  ^regions— border  ) } 

The  charges  and  currents  that  we  want  to  evaluate  in  region  II  can  be 
computed  using  the  corresponding  charge  and  current  functions  of  region  I: 


Qg(Vgs ,  Vbs ,  Vds)  =  Qg(Vgs,  5  Vds)', 


Qb(vgs,vbs,vds)  =  Qb(vgs,vbs,vdsy, 


Qd(Vgs,  14,  Vds)  =  Qs(Vgs ,  Vb3 ,  Vds)', 
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while 


Qs(vgs,vbs,vds)  =  QdIn(vga,vba,vdsy 


Ids(Vgs,  vba,  vds)  =  -Idt(Vga,  vba,  Vds). 

Special  care  has  to  be  taken  when  derivatives  are  calculated.  To  wit: 


9Qg 

dVds 


(V)  = 


9Qg 

dVn 


(V)  -  ^(V)  - 


gs 


dVbs 


ds 


The  procedures  evaluating  the  functions  were  written  as  the  above  im¬ 
plies:  Given  the  V’s,  the  V  is  computed  and  the  value  of  the  appropriate 
function  in  region  I  is  calculated;  the  signs,  etc.,  are  next  adjusted  according 
to  the  evaluated  function. 


3.3  Regions  III,  IV,  and  V 

We  know  next  to  nothing  of  the  behavior  of  the  intrinsic  transistor  outside 
regions  I  and  II;  nor  could  we  find  any  literature  on  it.  Since  we  needed  the 
functions  to  be  define  over  thewhole  domain,  we  extended  the  definitions  of 
the  basic  functions,  Ida,  Qs ,  Qb:  Qd,  and  Qg,  outside  regions  I  and  II  to  satisfy 
the  following  conditions:  the  functions  should  be  continuous  and  satisfy  the 
Lipschitz  conditions  over  the  entire  three  dimensional  space.  In  addition,  the 
extension  should  be  simple. 

Figure  6  describes  the  the  partition  of  the  space  into  regions  and  illus¬ 
trates  the  evaluation  of  Qd  and  Qs  in  each  region.  Given  a  point  {Vgs,  Vba,  Vds ) 
in  region  III,  IV,  or  V,  a  point  (Vgs,Vba,Vda)  is  calculated  on  the  boundary 
of  regions  I  or  region  II.  The  value  of  a  function,  one  of  Ida,  Qa,  Qb,  Qd,  or 
Qg ,  at  (Vga,  Vba,  Vds)  is  the  value  of  the  function  at  (Vga,  Vba,  Vda). 

The  definitions  of  the  regions  and  the  transformations  are  as  follows 

•  Region  III: 

^  ^regions— border  Vds  ^  0 ) ) 

Vg. =  v^s? 

—  ^regions— border  i 

vds  =  Vds. 


{(VSs,Hs,Vds)|(Hs 

Vba 
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•  Region  IV: 


{  ( hgs ?  Vbs  5  brfs  )  |  (  bf>s  ^  ^regions— border  )an(l( \’ds  ^  ^6 5  ^regions—  border  )  } 


I'gs  —  Vgs, 
bf>«  —  hfe*} 

V/s  —  iJ).s  ^regions— border  • 


•  Region  V: 


{  ( V/s  5  ^fcs  ?  Vis  )  |  (  V>s  ^  4>regions—border^&Tld(Vda  <C  0)} 


V/s  —  Igs5 


If)s  —  ^regions— border  i 

Vds  =  o 


Discussion:  No  claim  is  made  that  the  above  equations  capture  the  be¬ 
havior  of  the  intrinsic  transistor  in  regions  III,  IV  and  V.  As  will  be  shown  in 
the  sequel,  they  give  us  a  model  that  has  a  unique  solution  under  some  rea¬ 
sonable  conditions  and,  thus,  when  the  simulation  enters  these  region  we  can 
emit  a  warning  message  and  continue  instead  of  ending  with  a  programming 
bug. 

4  Graphical  presentation  of  the  results 

There  is  a  basic  difficulty  in  displaying  functions  of  three  variables  as  are 
the  current  and  charges  of  the  intrinsic  transistor.  To  gain  some  insight  into 
the  properties  of  these  functions  we  did  the  following:  Consider  figure  6  as 
a  plane,  Vgs  =  constant ,  in  the  space  (Vgs,Vbs,Vds).  Chose  a  value  for  Vb„ 
which  defines  the  dotted  line  (the  path)  in  the  figure.  Figures  10  to  18  depict 
the  current  Ids,  the  charges,  and  some  small  signal  parameters  commonly 
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used  by  analog  circuits  designers,  along  this  path.  The  region  boundaries  are 
shown  as  vertical  lines  in  the  figures. 

The  values  that  we  have  chosen  and  which  are  common  to  all  the  figures 
are  as  follows: 

Vbs  =  -5V; 

Vda  varies  continuously  from  -10V  to  10V; 

Vgs  is  a  parameter  value;  The  values  that  were  chosen  for  it  are  OV,  2V,  4V, 
6V  and  8V. 

The  following  parameters  were  chosen  for  the  representation  of  the  results: 
w  _  I 

L  1 

dox  =  700A  (Oxide  thickness) 

Na  =  3  •  1015  cm-3  (Doping  concentration) 

Lov  =  0.1 /i  (Channel  length  overlap) 

<f>ms  =  —0.93V  (Metal  semiconductor  potential) 

Qo  =  0  cm~3  (Effective  oxide  charge) 

Hon  =  1200  (Electron  mobility) 

Hop  =  250  ^  (Hole  mobility) 

Vsat  =  1  •  105  (Saturation  velocity) 

The  path  that  was  chosen  shows  the  behavior  of  the  functions  in  region  I  and 
II,  the  most  commonly  used  regions,  and  in  region  IV. 

5  Determining  the  model’s  characteristics  by  measure¬ 
ment 

One  major  property  that  we  expect  models  of  physical  devices  to  have  is  the 
ability  to  measure  experimentally  the  device  characteristics  from  the  device’s 
terminals,  and  then  to  use  the  characteristics  in  simulation.  The  SPICE  level 
3  model,  for  example  has  this  property.  In  this  section  we  present  a  way  of 
measuring  the  model  characteristics  from  the  terminals. 

Given  a  MOS  transistor  device,  we  would  like  to  extract  Qg,  Q &,  Qj,  Qs 
and  Ids  as  functions  of  Vgs ,  Vjs,  and  V&s.  Figure  4  shows  some  parasitic  capac¬ 
itors  in  addition  to  the  capacitors  of  the  intrinsic  model.  These  capacitors 
represent  the  capacitance  of  the  source-bulk  and  drain-bulk  pn  junctions, 
and  the  overlap  capacitance  between  the  gate  and  the  source,  and  between 
the  gate  and  the  drain. 

When  measuring  the  charge- volt  age  characteristics  for  each  of  the  tran- 
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sistor’s  terminal  pairs,  it  is  not  possible  to  distinguish  between  the  charge 
that  is  supplied  to  the  parasitic  capacitor  and  the  charge  that  is  supplied  to 
the  intrinsic  capacitor  that  appears  in  parallel  to  it  without  prior  assump¬ 
tions  on  the  behavior  of  the  voltage-charge  characteristics  of  each  one  of 
these  separately.  Thus,  in  this  paper  we  measure  the  two  charges  together. 
The  measurement  can  be  done  by  either  small  signal  capacitive  measure¬ 
ment  (C-V  measurement)  or  large  signal  charge  voltage  measurement.  In 
C-V  measurement  we  get  the  incremental  capacitance  and,  in  the  sequel,  it 
is  assumed  that  this  is  the  method  used. 

Let  Cij  be  the  total  incremental  capacitance  between  nodes  i  and  j  in  the 
the  device.  Thus,  we  get: 


CgS  —  Cpga 


+ 


dQg 


dV0 


gs 


vbs,vd, 


where  Cpgs  is  the  parasitic  capacitance  between  the  gate  and  the  source  (the 
overlapping  capacitance); 


Cgd  =  Cx 


pgd ? 


where  Cpgd  is  the  parasitic  capacitance  between  the  gate  and  the  drain  (the 
overlapping  capacitance); 


cds  = 


dQd 


dVda 


Vg„Vbt 


Cbs  —  Cpba 


dQb 

dVbs 


Cpbs  is  the  source-bulk  pn  junction  capacitance; 


Cbd  =  Cpdb 

Cpdb  is  the  drain  bulk  pn  junction  capacitance.  We  denote  by  CtJ  the  in¬ 
cremental  capacity  measurement  between  terminals  i  and  j  when  the  other 
terminals  are  connected  to  a  constant  DC  level. 

The  measurement  results  for  each  set  of  Vbs  Vgs  and  Vds  are: 

C gs  C gs  -|-  Cgd 

cbs  =  Cbs  +  Cu 
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Cds  —  Cds  +  Cdb  +  Cgd 
Cgd  ~  Cgd  +  Cdb  H"  Cds 
C gb  Cbs  4"  Cdb 

The  matrix  form  of  these  results  is: 


'1  1  0  0  0' 

'  C„  ' 

'  Cas' 

0  10  0  1 

Cgd 

Cbs 

0  1110 

Cds 

— 

Cds 

0  0  111 

Cdb 

Cgd 

0  0  0  1  1 

Cbs 

.  Cgb  . 

(4) 


We  would  like  to  obtain  the  incremental  capacitance  vector  vs.  the  measured 
capacitance  vector,  hence  we  get: 


’  Cgs  ' 

'  1 

-0.5 

-0.5 

0.5 

0 

Cgd 

0 

0.5 

0.5 

-0.5 

0 

Cds 

— 

0 

0 

0 

1 

-1 

Cdb 

0 

-0.5 

0.5 

-0.5 

1 

.  Cbs 

0 

0.5 

-0.5 

0.5 

0 

- 

' 

Cbs 

Cds 

Cgd 

- 

■Pi 

cr  ’ 

1 _ 

(5) 


Therefore,  we  can  get  the  incremental  capacitances  of  the  model’s  com¬ 
ponents  by  measurements.  In  order  to  get  the  characteristics,  charges  vs. 
voltages,  we  need  one  point  on  the  characteristics  in  addition  to  the  in¬ 
crements.  One  such  point  is:  For  Vgs  =  F/fc,  Vds  =  Vbs  =  0,  Qg  =  -Q0, 
Qd  =  Qb  =  Qs  =  0. 

These  measurements  can  be  used  to  compare  the  theoretical  derived 
model  with  the  actual  behavior  of  the  transistor,  or  to  obtain  the  charac¬ 
teristics  of  the  model  experimentally. 


6  Some  properties  of  the  of  the  model 

This  section  describes  some  qualitative  properties  of  the  model  all  of  which 
related  to  simulations.  These  properties  are  invertability,  satisfiability  of  Lip- 
schitz  conditions,  monotoneity  of  the  model  and  the  existence  and  uniqueness 
of  solution  of  networks  containing  such  models  as  components. 
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Is  the  model  charge  controlled?  We  are  accustomed  to  expect  the 
model  to  be  charge-controlled  [4].  The  expressions  defining  the  charges  and 
currents  in  our  model  are  all  expressed  in  terms  of  the  voltages  Vj,s,Vbs,and 
Vds ■  Does  the  function  which  defines  the  charges  in  terms  of  the  voltages 
have  an  inverse? 

Simple  inspection  of  the  intrinsic  model  yields  that  (a)  the  expressions 
in  regions  I  and  II  are  too  complicated  for  eye-ball  conclusions,  and  (b)  The 
functions  are  “flat”  on  the  other  regions  and,  therefore,  no  inverse  of  the 
intrinsic  model  can  exist. 

If,  however,  the  parasitics  are  added  to  the  picture  and  the  model  is 
considered  as  a  whole,  the  situation  is  not  obvious.  So,  we  consider  the  model 
with  its  parasitics,  i.e.,  the  linear  overlap  capacitors  Cpgs  and  Cpgd ,  and  thepn 
junction  capacitors  Cpbs  and  Cpbd •  Since  Cpbs  and  Cpbd  are  not  define  outside 
regions  I  and  II  they  were  extended  by  a  small  linear  capacitors  defined  over 
the  whole  domain.  Now,  the  answer  is  no  longer  obvious.  The  next  step  is 
the  use  of  a  global  inverse  theorem  and  an  experimental  verification  of  the 
conditions.  Let  us  choose  the  source  s  as  a  reference  node,  and  let  Qg,  Qb:am\ 
Qd  denote  the  node  charges  (cut-set  charges  associated  with  the  nodes)  which 
include  the  charge  of  the  parasitics.  Let  h  be  the  mapping  from  the  voltages 
to  the  charges,  i.e., 


The  model  is  charge-controlled  if  h  has  a  global  inverse.  To  prove  that  h  has 
the  desired  property  we  use  the  following  theorem. 

Palais’s  Theorem  [16]:  Let  h  be  a  differentiable  function  from  Rn  to  Rn. 
Then  h  is  a  diffeomorphism  3,  provided  that  h  satisfies  the  two  conditions 

lim 

INK00 

(it)  det  Dh(x )  /  0  for  all  x , 

were  Dh(x )  is  the  Jacobian  of  h  at  x. 

Our  functions  are  continuous  and  differentiable  inside  the  regions  but  are 
not  differentiable  at  the  regions  boundary.  On  one  hand,  this  implies  that 
if  an  inverse  exists,  it  would  not  be  differentiable;  On  the  other  hand,  (u) 

3h  has  a  continuous  differentiable  inverse  over  the  domain 


||/i(:r)||  =  oo 


19 


provides  local  homoeomorphism  and,  therefore,  we  have  to  prove  local  ho- 
moeomorphism  on  the  boundaries  some  other  way.  Now,  for  each  boundary, 
derivatives  on  the  “right”  and  the  “left”  exist.  Thus,  to  insure  local  ho¬ 
moeomorphism  it  is  sufficient  to  show  that,  at  each  point  on  the  boundary, 
the  corresponding  partial  derivatives  have  the  same  sign.  Only  the  normal 
derivatives  have  to  be  checked  as  continuity  implies  that  derivatives  along 
the  boundary  are  equal.  This  check  was  done  experimentally;  the  derivatives 
were  found  symbolically  and  evaluated  numerically  at  various  points  on  the 
boundaries. 

Let  V  be  the  vector  of  voltages  and  Q  -  the  vector  of  charges.  As  ||V|| 
sufficiently  increased,  so  does  ||C?||;  While  the  charges  of  the  the  intrinsic 
capacitors  do  not  increase  for  any  increase  in  V,  the  charges  of  the  linear 
capacitors  do  increase.  Thus,  condition  (?)  above  is  satisfied. 

Condition  (ii)  inside  the  regions  has  been  checked  numerically.  The  par¬ 
tial  derivatives  have  been  found  symbolically  and  the  determinant  was  cal¬ 
culated  at  various  points  in  space.  The  result  is  shown  by  figure  7. 

Thus,  strictly  speaking,  the  above  are  numerical  evidence  that  the  model 
is  charge-controlled  rather  than  mathematic  proofs.  A  proof  of  (ii)  can  be 
provided  by  the  evaluation  of  the  determinant  in  each  region  and  showing 
that  it  does  not  vanish. 

The  charges  and  Ijs  as  functions  of  the  voltages  satisfy  the  Lip- 
schitz  conditions  4. 

Since  all  the  above  functions  are  continuous  and,  on  each  region,  they  are 
defined  as  expressions,  the  above  is  determined  by  inspection. 

The  cutset  charges  are  strictly  monotonic  functions  of  the  volt¬ 
ages. 

A  function  /  from  Rn  to  Rn  is  called  monotonic  if  for  every  xi,x2  E 
Rn,x i  /  x2, 

(x!  -  x2,f(x x)  -  f(x2))  >  0. 

4 A  function  /  satisfy  the  Lipschitz  conditions  on  some  domain  D  if  there  exists  a 
k  >  0,  such  that  for  every  x\  and  x2  in  D 

||  /(*i)  -  f(x2)  ||<  k  ||  xi  —  x2  ||  . 

Satisfiability  of  Lipschitz  conditions  by  /  implies  existence  and  uniqueness  of  solution  of 
^ x  =  f(x).  For  differentiable  functions,  a  uniformly  bounded  derivatives  implies  that  the 
Lipschitz  conditions  are  satisfied. 
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THE  JACOBIAN  DETERMINANT 


Z  -  Jacobian  determinant 
X- Vds  (-10V- 10V) 

Y  -  Vgs  (0V- 10V) 


Figure  7:  The  jacobian  determinant 


In  the  above,  if  the  (>)  is  replaced  by  (>)  the  function  is  called  strictly 
monotonic.  If  the  right  handside  is  replaces  by  >  c||xi  —  a^H2  the  function  is 
called  strongly  monotonic. 

If  /  is  one-dimensional  and  differentiable,  the  above  means  that  the 
derivative  of  /  is  positive  if  /  is  strictly  monotonic  and  positive  or  zero 
if  it  is  monotonic. 

Our  interest  in  this  property  is  motivated  by  certain  fixed-point  theo¬ 
rems  on  monotonic  operators  [17,  18].  These  theorems  are  used  in  [15]  to 
derive  sufficient  conditions  for  the  existence  and  uniqueness  of  the  solutions 
of  nonlinear  networks  containing  transistors. 

The  characteristics  of  the  intrinsic  transistor  are  monotonic  but  not  strictly 
monotonic  since  their  their  derivatives  in  regions  II  IV  and  V  are  zero  (see 
figures  8  and  9).  Here  the  parasitic  capacitors  help  in  the  following  way: 
First,  the  overlap  capacitors  are  linear  and  satisfy  both  the  Lipschitz  condi¬ 
tions  and  the  monotonicity  conditions.  However,  the  value  of  the  pn  junction 
incremental  capacitors  is  zero  when  the  corresponding  diode  is  conducting. 
To  correct  that,  we  added  a  small  linear  capacitor  in  parallel  with  the  pn 
junction  capacitor  and  consider  the  sum  of  these  two  capacitors  as  the  par¬ 
asitic  capacitor.  Now,  if  the  transistor  is  considered  as  a  whole,  the  intrinsic 
model  plus  the  parasitics,  it  is  strictly  monotonic.  In  fact,  it  is  even  strongly 
monotonic  with  the  coefficient  c  determined  by  the  parasitics. 

7  Summary  and  Conclusions 

We  have  presented  a  model  which  is  defined  for  any  value  of  the  terminal 
voltages,  has  continuous  characteristics  (no  jumps)  all  of  which  satisfy  Lip¬ 
schitz  conditions.  The  model  is  defined  in  terms  of  physical  constants  and 
the  device  geometry.  Basically,  for  region  I,  the  model  summarizes  the  phys¬ 
ical  model  which  appears  in  [3].  Some  changes  were  introduced,  however,  to 
make  its  characteristics  continuous. 

The  model  is  summarized  as  a  non-linear  network  with  capacitors  rather 
than  with  incremental  capacitors  as  in  [4,  5,  6]  or,  in  fact,  SPICE  [7],  The 
translation  of  the  physical  model  equations  to  capacitors  is  natural,  although 
the  concept  requires  some  getting  used  to.  The  same  network  approach  can 
be  used  to  express  the  equations  of  the  short  channel  transistor  in  network 
terms.  Another  extension  which  is  conceptually  simple  is  the  modeling  for 
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A  V  *  A  Q 


Z- A  V*  AQ 
Y  -  Vgs  (OV  -  10V) 
X- Vds  (-10V-  10V) 
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higher  frequencies  by  dividing  the  transistors  to  cross-sections,  and  modeling 
each  such  section  with  a  lumped  network  model  similar  to  the  network  model 
presented  here. 

The  model  is  a  “large  signal  model”  and  it  suffices  for  simulation.  More¬ 
over,  conditions  for  existence  and  uniqueness  of  response  have  been  devel¬ 
oped. 

By  taking  appropriate  partial  derivatives  small  signal  behavior  can  be 
investigated  around  any  operating  point.  These  partial  derivatives  are  the 
incremental  capacitors,  incremental  resistors  and  gains.  The  partial  deriva¬ 
tives  are  continuous  everywhere  but  on  the  boundaries.  I.e.,  as  a  result  of 
continuity,  the  same  value  is  obtained  for  the  derivative  ‘along’  the  boundary 
when  the  derivative  is  evaluated  on  either  side  of  a  boundary.  The  normal 
derivative,  however,  evaluated  on  one  side  of  the  is  different  from  the  deriva¬ 
tive  evaluated  at  the  same  boundary  point,  but  which  is  considered  to  be  on 
the  other  side  of  a  boundary.  In  all  cases,  we  observed  that  the  two  deriva¬ 
tives  have  the  same  sign,  a  feature  used  in  the  proof  of  the  existence  of  the 
inverse. 

As  a  result  of  the  discontinuity  of  the  partial  derivatives,  the  incremental 
capacitors,  etc.,  are  discontinuous  functions  of  the  voltages.  The  technique 
used  in  [8]  can  be  used  to  smooth  the  derivatives  at  the  boundaries. 

It  seems  that  transistors  models  with  less  jumps  than  in  SPICE  [7]  are 
available  as  proprietary  models  in  industry  [9].  Since  these  models  are  private 
their  exact  properties  are  not  clear  to  us.  It  seems  that  [8]  and  ours  are  the 
only  “smooth”  models  publicly  available. 

Symbolic  analysis  played  a  major  role  in  our  work.  We  used  it  to  derive 
the  derivatives  of  the  charges  and  the  current  and  also  to  investigate  the 
properties  that  the  model  satisfies. 

The  model  parameters  can  be  measured  from  the  terminals.  Thus,  the 
model  provides  a  framework  for  experimentally  defined  models.  I.e.,  the 
topology  and,  maybe,  some  parameters  are  assumed  as  in  the  model,  and 
the  device  characteristics  can  be  measured  from  the  terminals  (see  [5]).  We 
have  not  conducted  such  measurements.  This  interesting  work  has  yet  to  be 
done.  Its  results  are  expected  to  shade  light  on  the  model’s  accuracy,  compare 
it  in  detail  with  other  models,  and  in  particular,  compare  the  multi-section 
network  model  with  a  transistor  operating  at  higher  frequencies. 
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8  Appendix  A:  Expressions  for  currents  and  charges 
in  region  I 

This  appendix  contains  the  expressions  for  the  long  channel  device.  Similar 
expressions  for  Ids  of  transistors  with  small  dimensions  can  be  found  in  [3]. 
We  did  not  find  in  the  literature  expressions  for  the  charges  of  small  dimen¬ 
sion  transistors.  We  have  looked  at  the  expressions  for  small  dimensions  Ids 
and  corrected  the  discontinuities  on  the  boundaries.  The  resulting  function 
obeys  the  Lipschitz  conditions  and  can  be  used  in  the  model  in  the  same  way 
as  the  long  channel  Ids- 

Some  common  physical  relations: 

Fermi  level  for  p-type  material:  4>f  =  ~ 


<f>s  =  2  (f>F 


C  — 

V-y  OX  - 


dn 


7  = 


es  q  Na 
— 

^  ox 


Vth  —  Vto  +  7  {^\J 4>s  —  vbs  —  y^Q 


where  Vtn  is: 


Vto  =  Vfb  +  4>s  +  7 
Ids  in  saturation  and  linear  regions: 

W 

Ids  =  Io  +  2"^ hoCoxVdsat{Vgs  —  —  «2)(  1  +  A  Vds) 


where  a  is  defined  as: 


Vu 


a  =  l--fV  for  Vds<Vdsat 

Vdsat 

Q  =  0  for  Vds  >  Vdsat 

Ids  for  the  sub-threshold  region,  i.e.,  Vgs  <  Vth- 

Vgs  -  Vth 


Ids  =  Io  exp 


3  kT 

2  q 


(6) 


(7) 
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where 


Io 


W 

T 


kT 

fi0Cox(—)21.5e 


-1 


qg  for  the  accumulation  region,  Vg8  <  (Vfb  +  Vbs): 


Qa  =  WLCox(Vgs  -  Vbs  -  Vfb)  -  Q0  (8) 


Qg  for  {Vfb  +  Vba)  <  Vgs  <  Vth: 

Qg  -  WLC°x-~j^\j4)s  -  Vbs  ~  Qo  (9) 

where  q  and  qvth  are: 


-  ~Y  +  \j^  +  Vgs  -  14s  -  Vfb 

qvth  =  ^  +  Vth  -  Vbs  -  Vfb 

Qg  for  the  linear  and  saturation  regions,  Vgs  >  Vth'- 

Q„  =  WLCo,  (vd„,(6  +  -  H,  )  -  <3o  (10) 

where 

<y=  7 

2\/  —  14s 

i/  _  fgs  — 

i  .  r 

1+0 

Qb  for  the  accumulation  region,  Vgs  <  {Vfb  +  14*): 

Qb  =  —WLCox{Vgs  -  Vba  -  V^)  (11) 

Qb  for  {Vfb  4-  14«)  <  14s  <  Vth- 

Qb  =  —WLCox  7 y/<f>.  -  Vbs  (12) 
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Qb  for  the  linear  and  saturation  regions,  Vgs  >  Vth: 

Qb  =  -  WLCox  [Vdsat  8(1  +  ^  Q  +  Q2)  + 

Qd- 

Qd  =  0  for  <  Vth- 
Qd  for  Vga  >  Vth: 

n  in/  t  n  dr  t  /  N  4  +  8a  +  12a2  +  6 a3 

Qi  =  —WLCox  (Vi.  -  V,t) - - 

Qs  =  0  for  14*  <  Vt/,. 

Qs  for  Vgs  >  Vth: 

n  in/  t r*  ii /  i/  \  6  +  12a  +  8a2  +  4a3 

Qs  =  -WLCox  (Vgs  -  Vth) - 7—  ,  ,2 - 

15(1  +  ay 

an  example  of  a  region  function  (for  Ids) 


•  Region  I 


•  Region  II 


where: 


Ids  =  f  (Vgs,  Vbs  ,  Vds) 


Ids=~f(Vgs,Vbs,Vds ) 


Vgs  =  Vgs  ~  Vds 
Vbs  =  Vbs  -  Vds 
Vds  =  -Vis 


Ids  =  f  (Vgs, 


regions— border 


Vds) 


Vt. 


•  Region  III 


•  Region  IV 


Ids  =  — /  (K*  >  Hs  ,  Vds ) 


where: 


Vgs  =  Vjs  "I"  ^regions— border  bfcs 

bf>s  —  ^regions— border 
Vds  —  ^regions— border 

•  Region  v 

Ids  —  f  (V/s  5  ^regions— border  i  0) 

An  example  of  a  derivative  function  (§yJL) 

for  Vgs  <  Vth 
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. 


9  Appendix  C:  Obtaining  the  Characteristics  Derivative- 
An  Example 

The  following  example  illustrates  how  the  derivatives  of  the  charges  have  been 
derived.  The  following  is  the  scheme  program  specifying  Qg  as  a  function 
of  the  voltages  and  the  physical  parameters.  It  is  followed  by  a  program 
that  calculates  the  partial  derivative  of  Qg  with  respect  to  Vgs  while  all  other 
parameters  are  constants.  The  resulting  program  is  the  third  item  in  this 
appendix. 
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(define  v->Qg 

(lambda  (  Vgs  Vbs  Vds  Vth  mu-zero  W  L  ) 

(let  ( (Basic-v->Qg 

(named-lambda  (Basic-v->Qg  Vgs  Vbs  Vds  Vth  mu-zero  W  L  ) 

(define  (fl  vgs  vbs  vds) 

(+  (*  W  L  Cox 

(let((q  (+  (-  (*  0.8  gamma)) 

(sqrt  (+ 

(*  gamma  gamma  0.25)  Vgs  (-  Vbs)  (-  Vfb))))) 
(qvth  (+(-(*  0.5  gamma)) 

(sqrt  (+ 

(*  gamma  gamma  0.25)  Vth  (-  Vbs)  (-  Vfb)))))) 
(/  (*  gamma  (sqrt  (-  Phis  Vbs))  q)  qvth))))) 

(define  (f2  vgs  vbs  vds) 

(let*  ((delta  (/  gamma  (*  2.  (sq-root  (+1.  (-  Phis  Vbs)))))) 

(Vdsat  (/  (-  Vgs  Vth)  (+  1.  delta))) 

(alpha  (if  (and  (<=  Vds  Vdsat)  (>  Vdsat  0.)) 

(-1.  (/  Vds  Vdsat));;  Linear  region 
0 .  ) ) ) ; ; saturation 

(*  W  L  Cox 
(+  (*  Vdsat 

(+  delta  (*  0.667  (/  (+1.  alpha  (*  alpha  alpha)) 

(+  1.  alpha))))) 

(*  gamma  (sq-root  (-  Phis  Vbs))))) 

)) 

(cond  ((<  Vgs  (+  Vfb  Vbs))  (+  (*  W  L  Cox  (+  Vgs  (-  Vbs)  (-  Vfb))))) 

((<  Vgs  Vth)  (fl  Vgs  vbs  vds)) 

(else  (f 2  vgs  vbs  vds))) 

))) 

(regions*  Basic-v->Qg  Vgs  Vbs  Vds  Vth  mu-zero  W  L)))) 

The  following  is  evaluated  in  the  symbolic  environment  to  yield  the  derivative  needed. 

A  special  utility  function  is  used:  ((pd  i)f)  finds  the  partial  derivative  of  a  scheme  func¬ 
tion  /  with  respect  to  the  i-th  variable.  Let  f  be  a  lambda  expression  of  two  variables, 

{{(pd  1  )/)'<*  ' b )  finds  the  function  with  respect  to  the  first  variable  and  then  finds  the 
value  of  the  derivative  at  the  (symbolic)  point  a, 6. 

(define  d/dvgs-Qg 

* (lambda  (  Vgs  Vbs  Vds  Vth  mu-zero  W  L  ) 

(let  ( (Basic-d/dvgs-Qg 
(lambda  (Vgs  Vbs  Vds  Vth  mu-zero  W  L  ) 


30 


(define  fl 

(lambda  (vgs  vbs  vds) 

(*  W  L  Cox 

(let*((q  (+  (-  (*  0.5  gamma)) 

(sqrt  (+ 

(*  gamma  gamma  0.25)  Vgs  (-  Vbs)  (-  Vfb))))) 

(qvth  (+  (-  (*  0.5  gamma)) 

(sqrt  (+ 

(*  gamma  gamma  0.25)  Vth  (-  Vbs)  (-  Vfb)))))) 

, (let*  ((q  (lambda  (vgs)  (+  (-  (*  0.5  gamma)) 

(sqrt  (+ 

(*  gamma  gamma  0.25)  Vgs  (-  Vbs)  (-  Vfb)))))) 

(qvth  'qvth) 

(exp  (lambda  (q) 

(/  (*  gamma  (sqrt  (-  Phis  Vbs))  q)  qvth)))) 

(*  ( C (pd  1)  exp)  'q)  (((pd  1)  q)  'vgs))))))) 

;  cutoff  270 

(define  f2 

(lambda  (vgs  vbs  vds) 

(let*  ((delta  (/  gamma  (*  2.  (sqrt  (-  Phis  Vbs))))) 
(Vdsat  (/  (-  Vgs  Vth)  (+  1.  delta))) 

(alpha  (if  (and  (<=  Vds  Vdsat)  (>  Vdsat  0.)) 

(-  1.  (/  Vds  Vdsat));;  Linear  region 
0 .  ) ) ) ; ; saturat ion 

(*  W  L  Cox 
(if  (>  Vdsat  0.) 

, (let*  ((delta  'delta) 

(vdsat  (lambda(Vgs)  (/  (-  Vgs  Vth)  (+  1  delta)))) 
(vdsat 1  '(/  (-  Vgs  Vth)  (+  1  delta))) 

(alpha  (lambda (vds at)  (-  1  (/  Vds  Vdsat)))) 

(exp  (lambda  (vdsat  alpha) 

(+  delta 
(*  2/3 

(/  (+  1  alpha  (*  alpha  alpha)) 

(+  1  alpha)))))) 

(expl 
'  (+  delta 
(*  2/3 

(/  (+  1  alpha  (*  alpha  alpha)) 
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(+  1  alpha)))))) 

(+  (*  expl  ( ( (pd  1)  vdsat)  'vgs)) 

(*  vdsat 1  (((pd  2)  exp)  ' vdsat  'alpha  ) 

(((pd  1)  alpha)  'vdsat) 

(((pd  1)  vdsat)  ' vgs ) ) ) ) 

(/  1.  (+1.  delta))) 

)))) 

(cond  ((<  Vgs  (+  Vfb  Vbs))  (*  W  L  Cox  ))  ; Accumulation 
((<  Vgs  Vth)  (fl  Vgs  vbs  vds)) 

(else  (f 2  vgs  vbs  vds)))) 

)) 


(regions*  Basic-d/dvgs-Qg  Vgs  Vbs  Vds  Vth  mu-zero  W  L)))) 

The  result; 

(define  d/dvgs-qg 
(lambda 

(vgs  vbs  vds  vth  mu-zero  w  1) 

(let 

( (basic-d/dvgs-qg 
(lambda 

(vgs  vbs  vds  vth  mu-zero  w  1) 

(define 
f  1 

(lambda 
(vgs  vbs  vds) 

(* 

w 

1 

cox 

(let* 

CCq 

(+ 

(-  (*  .5  gamma)) 

(sqrt  (+  (*  gamma  gamma  .25)  vgs  (-  vbs)  (-  vfb))))) 
(qvth 
(+ 

(-  (*  .5  gamma)) 

(sqrt  (+  (*  gamma  gamma  .25)  vth  (-  vbs)  (-  vfb)))))) 

(* 

(/  (*  gamma  (sqrt  (+  phis  (*  -1  vbs))))  qvth) 

(/ 
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1/2 

(sqrt  (+  (*  .25  (expt  gamma  2))  (*  -1  vbs)  (*  -1  vfb)  vgs)))))))) 
(define 
f  2 

(lambda 
(vgs  vbs  vds) 

(let* 

((delta  (/  gamma  (*  2.  (sqrt  (-  phis  vbs))))) 

(vdsat  (/  (-  vgs  vth)  (+  1.  delta))) 

(alpha 

(if  (and  (<=  vds  vdsat)  (>  vdsat  0.))  (-  1.  (/  vds  vdsat))  0.))) 

(* 

w 

1 

cox 

(if 

(>  vdsat  0. ) 

(+ 

(* 

(+  delta  (*  2/3  (/  (+  1  alpha  (*  alpha  alpha))  (+  1  alpha)))) 

(/  1  (+  delta  1))) 

(* 

(/  (-  vgs  vth)  (+  1  delta)) 

(* 

(* 

(/ 

(  +  (*  2/3  (expt  alpha  2))  (*  4/3  alpha)) 

(+  (expt  alpha  2)  (*  2  alpha)  1)) 

(/  vds  (expt  vdsat  2))) 

(/  1  (+  delta  1))))) 

(/  1.  (+  1.  delta))))))) 

(cond  ((<  vgs  (+  vfb  vbs))  (*  w  1  cox)) 

((<  vgs  vth)  (fl  vgs  vbs  vds))  (else  (f2  vgs  vbs  vds)))))) 

(regions+  basic-d/dvgs-qg  vgs  vbs  vds  vth  mu-zero  w  1)))) 
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Figure  13:  Qd  for  Vbs=-5V 


gmb  vs.  Vds 


Figure  16:  gmb  for  Vbs=-5V 

gds  vs.  Vds 


Vds 

Figure  17:  gds  for  Vbs=-5V 
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d/dvds-Qg 


Figure  1$:  for  Vbs=-5V 
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